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ABSTRACT 



During the evolution of first stars, the CNO elements may emerge on their surfaces due to the mixing processes. 
Consequently, these stars may have winds driven purely by CNO elements. We study the properties of such stellar 
winds and discuss their influence on the surrounding environment. For this purpose, we used our own NLTE models 
and tested which stellar parameters of the first stars at different evolutionary stages result in CNO winds. If such winds 
are possible, we calculate their hydrodynamic structure and predict their parameters. We show that, while the studied 
stars do not have any wind driven purely by hydrogen and helium, CNO driven winds exist in more luminous stars. On 
the other hand, for very hot stars, CNO elements are too ionized to drive a wind. In most cases the derived mass-loss 
rate is much less than calculated with solar mixture of elements. This is because wind mass-loss rate in present hot 
stars is dominated by elements heavier than CNO. We conclude that, until a sufficient amount of these elements is 
created, the influence of line-driven winds is relatively small on the evolution of hot stars (which are not close to the 
Eddington limit). 
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1. Introduction 

The only elements created in non negligible amounts during 
primordia l nucle osynthesis are helium and hydrogen (e.g., 
ICoc et all I2004D . Consequently, first stars, which formed 
from the pristine gas processed during this primordial nu- 
cleosynthesis, were pure hydrogen-helium stars. Numerical 
simulations of first star formation show that the primordial 
nucleosynthesis not only determined the chemical compo- 
sition of these stars, but the absence of heavier elements 
also influenced the initial mass function of these st ars (e.g., 
iBromm et allll999t iNakamura k Umemural I2002D . In the 
absence of heavier elements, the only efficient cooling pro- 
cesses during the collapse of primordial haloes are those 
of molecular hydrogen. Since the molecular hydrogen line 
cooling is less efficient than, say cooling from dust, the tem- 
perature in the first star-forming regions was much higher 
than in the present-day star-forming regions. In the absence 
of fragmentation this led to very high masses of firs t stars 
(of the order of 10 - 100 M^. lOi^ukai fc Pallal[200.1^ . 

The general picture of the evolution of first stars is 
in many a spects different from the evolut i on of present 
hot s tars (jMarigo et al.l 120011 : iKlapp et all 120051 : iHirschil 
l2007t ). One of the most important differences is the pos- 
sible existence of pair-instability supern ovae, which create 
a ver y typical pattern of nuclear yields (jHeger fc WooslevI 
I2OO2I ). The possibility of creating a pair- instability super- 
nova is closely connected to the total amount of mass lost 
by an individual star during its evolution. If the star does 
lose a significant fraction of the mass during its evolu- 
tion, it may avoid th e pair-instability supernova eruption 
(jEkstrom et al.l 120071 ). 

Curren t massive stars lose their mass via line-driven 
winds (see lKudritzki fc Pulsl[2000l : lKrticka fc Kubatll2007aL 



for reviews). These winds are accelerated mainly due to the 
absorption in the lines of heavier elements, such as iron, 
carbon, nitrogen, and oxygen. For these stars the amount 
of mass lost by the star per unit of time (the mass-loss 
rate) can be in principle derived both from observations and 
from theoretical modelling (e.g. JPuls et al. 2006). However, 
given that there have not been any available observations 
of first stars up to now, we have to rely on the theoretical 
predictions alone in the case of first stars. 

Stellar winds of pure hydroge n-helium extremel y mas- 
sive first stars we re studied by iKrticka fc Kubafj (|2006l 
hereafter iPaper 13 ). We showed that the homogeneous line- 
driven winds of these stars are unlikely and that only an 
extremely weak pure hydrogen wind could be possible for 
stars very close to the Eddington limit. Stellar wi nds of very 
low-metallicity massive stars were studied by iKudritzkH 
(2002,) assuming the solar mixt u re of e lements. On the basis 
of these models, iMarigo et al.l (|2003D concludes that line- 
driven winds influence the evolution of non-rotating, very 
low-metallicity stars only negligibly. Stellar winds of stars 
at extremely low metallicity might be po ssible only if the 
star is very close to the Eddington limit (iKudritzkil l200d : 
IVink fc de Koteill2005HGrafener fc Hamannll2008l) 



During the core He-burning phase of fi rst stars, pri- 
mary nitrogen emerges on the stellar surface (|Mevnet et al.l 
2006; Hirschi 2007). This primary nitrogen is synthetised in 
the H-burning shell due to the rotational mixing of carbon 
and oxygen produced in the helium core. M oreover, evolu- 
tiona r y mode l s of massive first stars (e.g., iMevnet et al.l 
I2OO6I : iHirschil I2007D show that subsequent generation of 
stars could contain relatively large amounts of CNO ele- 
ments, whereas the abundance of iron-peak elements was 
very low. This picture is also supported by observations of 
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Table 1. Atoms and ions included in the NLTE calcula- 



tions. 


Ion 


Levels^^ 


Data 


Comment'^ 


Hi 


9 










Hii 


1 










He I 


14 






< 


70 000 K 


Hen 


14 










Hem 


1 










Ci 


26 




TcB 


< 


30 000 K 


Cii 


14 




Tcff 


< 


35 000 K 


Cm 


23 


OSTAR2002^ 


Tcs 


< 


70 000 K 


Civ 


25 


OSTAR2002 








Cv 


11 










Cvi 


10 




TcB 


> 


30 000 K 


Cvii 


1 




TcB 


> 


50 000 K 


Ni 


21 




TcB 


< 


30 000 K 


Nil 


14 




TcB 


< 


35 000 K 


Niii 


32 


OSTAR2002 


TcB 


< 


70 000 K 


Niv 


23 


OSTAR2002 








Nv 


16 


OSTAR2002 








Nvi 


15 




TeS 


> 


20 000 K 


Nvii 


1 




TeB 


> 


30 000 K 


Oi 


12 




TcB 


< 


30 000 K 


On 


50 




TcB 


< 


35 000 K 


Oiii 


29 


OSTAR2002 


TcB 


< 


70 000 K 


Oiv 


39 


OSTAR2002 








Ov 


14 










Ovi 


20 


OSTAR2002 


TcB 


> 


20 000 K 


Ovii 


1 




TcB 


> 


30 000 K 



^ An individual level or a set of levels merged into a superlevel. 
^ Note: the comment (if any) shows for which stars the individual 
ion is taken into an account. If the comment is missing, then the 

ion is considered for all modelS; 

^ The ionic model taken from lLanz fc Hubenvl lj2003l ). 



carbo n-rich, extremely metal-poor stars (e.g.. iNorris at all 
I1997D . 

However, the detailed study of stellar winds driven 
purely by CNO elements is still missing. Here we study 
such winds for the parameters typical of hot first stars. 

2. Description of models 

2.1. Basic assumptions 

The mod els used in this paper are based on the NLTE wind 
models of iKrticka fc Kuball (|2004 hereafter [Paper t ). Here 
we only summarise their basic features and describe im- 
provements in atomic data. 

We assume spherically symmetric stationary stellar 
wind. The excitation and ionization state of the considered 
elements is derived from the statistical equilibrium (NLTE) 
equations. For some ions we improved the atomic data 
and adopted the ionic model s from the OSTAR2002 grid 
of model stellar atmospheres (|Lanz fc Hubenvl [20031 . l2007l 
see Table [T]) . These io nic mod els are based mainly on the 
Opac i ty Project data j Peach e t al. 198 8; L uo fc Pradhan 
1989t iTullv et al.l Il990t ISeaton et all Il992l: iFernlev et all 
1999() . For some of these stars, the ionization fractions of 
some ions in the envelopes are so low that it is not meaning- 
ful to include these ions in the solution of the NLTE equa- 
tions. The conditions for including such ions are also given 
in Table [H The ionic list is modified in such a way that the 
highest considered ion is only included as its ground state. 



Table 2. Radius _R», mass M, and the effective temperature 
Teff of studied model stars 





ZAMS Stars 






Evolved stars 




Model 




M 


Tcff 


Model 




M 


Tcff 




[R©] 


[Mo] 


[kK] 




[R©1 


[Mq] 


[kK] 


M999 


4.23 


100 


94.4 


M999-0 


8.2 


100 


73.6 


M700 


3.44 


70 


89.5 


M999-1 


56.4 


100 


29.9 


M500 


2.82 


50 


84.1 


M999-2 


125 


100 


20.1 


M300 


2.10 


30 


74.0 


M999-3 


510 


100 


10.0 


M200 


1.65 


20 


65.3 


M500-1 


11.1 


50 


50.0 


M150 


1.48 


15 


57.3 


M500-2 


33.7 


50 


29.9 


M120 


1.42 


12 


49.9 


M500-3 


72.0 


50 


20.6 


MlOO 


1.37 


10 


44.5 


M500-4 


303 


50 


10.1 


M090 


1.34 


9 


41.6 


M200-1 


4.1 


20 


50.0 


M080 


1.31 


8 


38.5 


M200-2 


19.9 


20 


24.5 


M070 


1.30 


7 


34.8 


MlOO-1 


11.4 


10 


20.2 


M060 


1.27 


6 


31.4 


MlOO-2 


45.6 


10 


9.8 


M050 


1.23 


5 


27.7 


M050-1 


5.1 


5 


20.1 


M040 


1.17 


4 


23.6 










M030 


1.11 


3 


19.1 










M020 


1.00 


2 


13.7 











As in [Paper j . the radiative transfer problem is arti- 
ficially split into two parts, namely the radiative transfer 
in continuum and the radiative transfer in lines. The so- 
lution of the radiative transfer equation in continuum is 
based on the Feautrier method in the sp herical coordinates 
([Mihalas fc Hummer! [1971 iKuball [19931 ) with inclusion of 
all free-free and bound-free transitions of model ions, how- 
ever neglecting line transitions. 

The radiati ve tran sfer in lines is solved in the Sobolev 
approximation ([Castor 1974) neglecting continuum opacity 
and line overlaps. The radiative force is calculated in the 
Sobolev appr oximatio n using data from th e VALD database 
([Piskunov et al. 1995: i Kupka et al.lll999l) . To test the com- 
pleteness of our line list, we compared t hese data with the 
line-lists pro vided by t he NI ST database ([Wiese et al.ll996f) 
and given by iKurucd ([1992f ) and enlarged the original data 
set when necessary. The radiative force due to the light 
scattering on free electrons is also included in the model 
calculations. 

The surface emergent flux (i.e., the lower boundary con- 
dition for the radiative transfer in wind) is taken from H-He 
spherically symmetric NLTE model stellar atmospheres of 
,Kubat (2003, and references therein). 

Parameters of these stars (given in Table ^ were ob- 
tained according to the evolutionar y calculation s of ini tially 
zcro-metallicity stars derived by Marigo "eTal] ([200lh . We 
selected stellar parameters corresponding to the zero-age 
main sequence of these models and to later evolutionary 
phases to cover a larger area of the HR diagram. 

2.2. The wind tests 

Since we do not know in advance whether the radiative force 
is strong enough for a given star and chemical composition 
to drive a wind, we first tes t whether a wind may even ex- 
ist. For these tests (see also [Paper Ilf ). the hydrodynamical 
variables (velocity, temperature, and the density) are kept 
fixed. This enables us to calculate the model occupation 
numbers and the radiative force even in the case when the 
wind does not exist. The comparison of calculated radiative 
force and the gravitational force serves as a test of whether 
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the wind exists; i.e., the winds are only possible when the 
magnitude of the radiative acceleration g'^^'^ is greater than 
the magnitude of the gravity acceleration 5, 



(1) 



Besides the cotribution of lines, we also include the radia- 
tive acceleration due to the bound-free and free-free tran- 
sitions in the calculation of the total radiative acceleration 

^rad 

The velocity structure of our models is given b y an ar- 
tificial velocity law (see the discussion in iPaper Ij ) 



5 kTcs 
3 toh 



-I- 2 X 10^ cms" 



(2) 



where T^g is the stellar effective temperature, i?* the stel- 
lar radius, and mu the mass of the hydrogen atom. The 
density structure is obtained from the equation of continu- 
ity. In these models we assume a constant wind temper- 
ature O.SToff, and the electron density is consistently cal- 
culated from the ionization balance. Since for these wind 
tests we do not solve the equation of motion, it is nec- 
essary to specify wind mass-loss rate for which the wind 
existence is tested. For each set of stellar parameters the 
wind existence is tested for mass-loss rates 10~^ M0 year~^, 
lO~i°M0year-\ and IQ-^^ Mq year'^ 

To obtain a correct surface flux we accounted for the 
Doppler effect during the calculation of radiative transfer 
in lines, i.e. we shifted th e stellar sur face flux according to 
the actual wind velocity (|Babellll99^ . 

2.3. Hydrodynamic wind models 

For those stars, for which the wind test showed that the 
wind is possible, we also cal culated hydrodynamic wind 
models (described in iPaperT ). In these models we consis- 
tently solved the hydrodynamical equations, i.e. the con- 
tinuity equation, the momentum equation, and the energy 
equation with radiative cooling and heating included using 
the electrons thermal balance method (Kiibat e t al. 199!^. 

These models enable us to derive the radial variations 
of wind density, velocity, temperature, and the occupation 
numbers of individual excitation states. They especially en- 
able us to predict the wind mass-loss rate and the terminal 
velocity. 



these elements is by three to four orders of magnitude lower 
than the absolute value of the gravitational force. The most 
significant contribution to the radiative force is caused by 
the light scattering on free electrons for these stars. 

For cooler stars, helium may be partly singly ionized; 
however, these stars have very broad photospheric lines, 
making the flux at the line positions so low that it does not 
allow the wind to accelerate. For stars with T^s < 20 000 K, 
helium may in addition become partly neutral (in the en- 
velopes with highest densities). However, the flux in the 
ultraviolet part of the spectrum where the neutral helium 
resonance lines appear is so low that the radiative force due 
to these lines is too low to accelerate the wind. In all these 
envelopes the hydrogen is nearly completely ionized, so that 
the radiative force due to hydrogen lines is negligible. 

4. Analytical models of CNO line-driven winds 

Before discussing more detailed numerical NLTE wind 
models of CNO line-driven winds, we first pro vide sim plified 
analytic models of these wirids (c. f., lAbbott 1980: O wockil 
[200I iMdl^eier fc Shlos^[200l . 

The radiative acceleration in the Sobolev approxi mation 
is the sum of the contributions by individual lines (|Castoij 
I197I 

^rad = ^^Y, '^rjH,{,y,,) f ^i {\ + a/i^) (l - e"^") , 

^ linos 

(3) 

where Hc{vij) is the emergent flux from the stellar atmo- 
sphere at the frequency Vij of a given line, v is the radial 
wind velocity, p the wind density, 

fx— ^-1, (4) 



V dr 



.2\ 1/2 



(5) 



and the Sobolev optical d epth is given by (|Castoilll974l : 
iRvbicki k Hummeil[T978l) 



7re 



(6) 



mcVij \ gr gj J u (H- up^) ' 

where n^, rij, gi^ gj are the number densities and the statis- 
tical weights of individual states giving rise to a given line 
with the oscillator strength j . 



3. Pure H-He winds 



In iPaper Ij we showed that the extremely massive (Af > 
100 Mq), hot, pure hydrogen-helium stars do not have any 
wind. Stars very close to the Eddington limit with F > 
0.859 might be the only exception to this, since these stars 
can have very weak pure hydrogen wind. The absence of 
winds is connected with the material in the vicinity of these 
stars beeing ionized and consequently the line transitions 
are not strong enough to drive a wind. On the other hand, 
the material in the envelopes of the cooler, less massive stars 
studied here may be less ionized, thus giving the possibility 
to drive a wind. Here we test whether studied stars may 
have winds driven purely by hydrogen and helium. 

For very hot stars, with Tcff > 5 • 10* K both hydrogen 
and helium are completely ionized in the envelopes; con- 
sequently, the radiative force due to the line transitions of 



4.1. Wind driven by optically thick lines 

Generally, the stellar wind of present hot stars is driven 
by an ensemble of optically thick and thin lines, the op- 
tically thick lines being important due to their strength 
and the optically thin ones due to their large number (e.g., 
iPuls et al]|2000l ). The situation in the winds driven purely 
by CNO elements is different. Contrary to say, iron, these 
elements do not have many lines to drive wind, and the op- 
tically thick lines become much more important. In many 
cases there are just few optically thick lines that drive a 
wind. In such cases (for ^ 3 to obtain a precision of 
5%), the radiative acceleration Eq. ^ is given roughly by 

V dr 2 V 1; dr / 



rad 
ythick 



,2 



E 



(7) 
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The momentum equation neglecting the gas pressure term 
and introducing the Eddington parameter, 



rad 



r = 



1 CTeL 



(8) 



5thfck 



(1 - T)GM 



(9) 



takes the form 

dv dv 
F{r,v, — ) = V— 

dr dr i- 

Here, gl'^'^ is the radiative acceleration due to the free 
electrons, L the stellar luminosity, CTo = naSc/p, Sc the 
Thomson scattering cross-section, Ue the electron number 
density, and g = GM/r"^ . The momentum equation Q has 
the critical point defined by dF{r,v, ^)/d ^ = 

(rcrit is the critical point radius) 



M = Airr pv = 



1 - 



thick 
lines 



VijHc{vij), 



which yields the wind mass- loss rate. Con t rary t o the orig- 
inal analysis presented by ICastor et al.l (|l975l . hereafter 
CAK), the critical point condition is simpler, as the radia- 
tive force depends on the velocity gradient linearly here. It 
would be possible to use the regularity condition to derive 
Tcrit and M (as was done in CAK); however, as in most 
cases Tcrit ~ -R*, the wind mass- loss rate can be roughly 
given by 



M 



(11) 



This equation states that the mass-loss rate of the wind 
driven purely by optically thick lines is approximately given 
by the photon ma ss- loss rate (L/c^) multi plied by the num- 
ber of thick hues (|Lucv fc Solomonlll970tlOwocki|[l99l . It 
can be shown that Eq. (jlOp also defines the point where 
the wii id velocity is equal to the speed of s o-called Abbott 
waves (|Abbottlll98(l : iFeldmeier et al.ll2008h . 

Using the mass-loss rate estimate Eq. ([Tl]) and the defi- 
nition of the mass- loss rate (M = Airr^pv), the momentum 
equation Eq. ([9]) can be rewritten as 



^_g^d.^(i-r)GM_gy^ (12) 

d?' 



This equation has a critical point at Tcrit — R*- To ensure 
the finiteness of the velocity derivative at this point, the 
righthand side of Eq. (fT^ should be equal to zero. This 
yields the value of the wind velocity at the critical point 



l^crit 



'(l-r)GM 



(13) 



The application of the I'Hopital rul e to the momentum 
equat ion Eq. at the critical point (|Lamers fc Cassinellil 
I1999D gives the value of the velocity derivative at this point: 



dv 



dr 4i?» 



1 /(i-r)GA/ 



(14) 



4.2. Theoretical limits for CNO line-driven winds 

In an accelerating outflow, the radiative force and the ther- 
mal pressure have to overcome gravity and inertia at each 
point in wind. As the thermal pressure is only effective at 
the base of the wind, the necessary condition for launching 
line-driven stellar wind is that the magnitude of radiative 
force should exceed the magnitude of gravitational force at 
a certain point in the wind, Eq. ([1]). We apply this condi- 
tion to derive the minimal mass-loss rate of homogeneous 
wind (i.e. composed of all constituents, namely hydrogen, 
helium, carbon, nitrogen, and oxygen). 

For a given set of lines, the radiative force is maximum if 
the lines are not self-shadowed, i.e. if the lines are optically 
thin and the line optical depths t < 1 (Ga vlev 1995,) . Thus, 
for given occupation numbers, t he line ac celeration Eq. ([3]) 
cannot be higher than (see also iPaper III . Eq. (2)) 



^rad, 



prricC^ 



' ^ lines 3] J 



(15) 

The condition of th e wind existence Eq. ([T|) then reads as 
(c.f. lKudritzkil[2003 ) 



_rad, 



9 



+ 9T={Q + l)9.>9 



(16) 



where the par ameter Q = g ' max^^rad introduced by 
iGavlevI ()1995f ). Neglecting the term with the number den- 
sity of the upper level in Eq. (|15p (which means neglect- 
ing stimulated emission) and substituting g ~ GM/r'^, the 
wind condition (1161) can be rewritten as 



lines 



-r > 1. 



(17) 



For a very low density of CNO elements, only the 
strongest lines (i.e., the resonance lines of the most abun- 
dant ions) significantly c ontribute to the radiative force 
(e.g., iKrticka et al.l [20061 ). In such a case, condition (fT7|) 
can be rewritten as 



r> 1, 



(18) 



where Zd is the mass fraction of element corresponding to 
the transition i ^ j (the ratio of the element density to the 
wind density), TOqi is the element mass, and the summation 
goes over the resonance lines of only the most abundant 
ions. If there is only one such element, then the minimum 
mass fraction of this element needed to drive the wind is 
given by 



In scaled quantities, Eq. (fT9|) takes the form of 



(19) 



0.04Aei (1 - F) 



E 



M \ f R* 
IM^J [Tr^ 



erg cm" 



(20) 
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where ^ci is the elemental mass in units of hydrogen mass. 
Consequently, the minimum metallicity necessary for driv- 
ing a wind depends on the surface gravity g = GM/R^, 
and on the stellar effective temperature that determines 
both the ionization balance and the radiation flux. 

The minimum mass-loss rate that can be driven from 
the star is given by the condition that at least one line 
remains to be optically thick at the critical point, i.e. t^j > 
1, where is given by Eq. ([6|). Using Eqs. (fT3l) and ^4]) . we 
can show (assuming rcrit ~ R* and for radial rays = 1) 
that the mass-loss rate should fulfil the condition 

TOelTTT-o — T)GM 



M > 



(21) 



6' fijZol 

For a typical case where there is only one optically thick 
line left, it can be readily shown (using the mass-loss rate 
estimate Eq. (fTTj) ) that the condition Eq. ((2T|) 

^^'''iy!' gc(M,,)/.,+r>i, (22) 

is nearly the same as the condition Eq. written for one 
line. Consequently, if there is a line capable of overcoming 
the gravity (according to Eq. (fTS)) '). it is automatically op- 
tically thick one, and this line sets the mass-loss rate (and 
vice versa). In other words, for very low-density winds the 
condition of maximum radiative force Eq. (|16p and the crit- 
ical point analysis Sect. 14.11 give the same results for the 
existence of the wind and the wind's mass-loss rate. 

5. CNO-driven wind models 

For the parameters of model stars given in Table [21 we 
tested whether these stars could have CNO driven winds. 
We tested the possibility of the wind for the mass-loss rates 
10-12 M0year-i, 10"^° Mq year"!, and 10"^ Mq year"!, 
and for the mass fraction of CNO elements in the range 
of 10"'' — 0.1 with a relative step of 10. We assumed the 
same number density for carbon, nitrogen, and oxygen. 

The regions in the HR diagram corresponding to the 
parameters of stars for which the wind tests are positive 
are given in Fig. [TJ The minimum mass fraction of CNO 
needed to drive the wind is also indicated there. 

For parameters for which the wind may exist, we calcu- 
lated hydrodynamic wind models and predicted the mass- 
loss rate (again assuming the same number density of 
CNO). Calculated wind parameters for individual stars are 
given in Appendix in Table ETT] The mass-loss rate of stud- 
ied stars can be fitted by the formula 

a2(logg-l-a3 log-L) 
M(L,Z,Teff) =ao 10log2(a4-ha5logL) + l+«'='°S^<=«^ 

(23) 

where the luminosity L is expressed in the solar luminosity 
units Lq, Z is the mass fraction of heavier elements, and 
Tgff is expressed in the units of Kelvin. The values of pa- 
rameters ap, . . . , ae fitted using the subroutine VARPRO 
(iGolub &: Perevrall97alBolstadll977D for the stellar groups 
with different effective temperatures are given in Table [31 
For log Z{a4 + as log L) + 1 < the predicted mass-loss 
rate is set to zero. 

In present-day hot stars, the modified wind momentum 

1/2 

Mvqc (i?/R(;)) d epends mainly on the stellar luminosity 
(jKudritzki fc Pul"il2000l and references therein) . According 
to its value, these stars can be divided into four groups (see 
Fig. [ID. 
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Fig. 1. The HR diagram with minimal metallicity needed 
to drive a wind. Stars for which the winds are possible (for 
a studied range of metallicities) are denoted by crosses, and 
stars for which the winds are not possible are denoted by 
circles. The dashed lines denote boundaries separating stars 
with a different minimum metallicity necessary for driving 
a wind. The boundaries were derived using Eq. (|18p and 
the Wien law . Ove rplotted are the evolutionary tracks of 
iMarigo et all (|200lD . 
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Fig. 2. The modified wind momentum-luminosity relation- 
ship for stars with different effective temperatures and 
Z = 0.01. Stars with TeS > 70 000 K do not have any wind. 
Overplotted is the mean dependence of the modified wind 
momentum predicted for O stars with solar metallicity by 
.Krticka fc Kubat (2007b,). 



5.1. Stars with no wind: Tefr> 70 000 K 

For the hottest stars with the effective temperature Toff > 
70 000 K, we did not find that the wind exists. Although 
these stars are very luminous, their envelope is highly ion- 
ized with ions VI-VII being the dominant ionization states 
of CNO elements. Since the factor (M/1 M©) (i?,/l R©)"^ 
in Eq. (j^O)) is larger than one for all these stars, a strong 
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Table 3. The best- fit parameters of Eq. for individual stellar groups 



Tcff [K] 


ao [Moyear"-^] 


Ql 


a2 




Q4 


as 




< 28 000 


1.19 X IQ-^^ 


0.8437 


0.2936 


0.7167 


0.1208 


0.0064 


2.4369 


28 000 - 40 000 


8.92 X IQ-" 


0.7895 


0.6123 


0.8181 


0.2848 


-0.020158 


-0.7394 


> 40 000 


2.39 X 10" 


0.5722 


0.8348 


0.7995 


0.3983 


-0.03542 


-5.363 



flux at the line position is needed to accelerate the wind. 
However, ions vi- vii have strong resonance lines in the far 
UV region (for frequencies higher than the He ii ionization 
threshold), where very low flux is emitted, so we do not find 
any wind for these stars. 



5.2. Hottest stars with winds: 40 000 K < Teff < 70 000 K 

Stars with the effective temperatures in the range 
40 000 K < Teff < 70 000 K may have a line-driven wind ac- 
celerated mainly by the Ov and O VI lines. The hydrogen- 
like ion Cvi and the helium- like one Nvi, which arc dom- 
inant ionization fractions of carbon and nitrogen, do not 
contribute significantly to the radiative force. The latter 
ions do not have enough strong resonance lines at the wave- 
lengths at which strong flux is emitted. The minimum mass 
fraction of CNO elements needed to drive a wind agrees 
with analytical estimate Eq. (PO)) and is relatively high, of 
the order of 10""^ for ZAMS stars and as low as 10"'* for 
more evolved stars with larger radii. 

The winds of hottest stars in this sample are accel- 
erated mainly by the Ovi resonance doublet 2s ^S® — > 
2p^P°. The analytic formula Eq. ([TT|) gives the asymp- 
totic value of the mass-loss rate if both lines of the dou- 
blet are optically thick (i.e., for high metallicity) , while 
for lower metallicity the lines become optically thinner and 
the mass-loss rate decreases. For example, for the model 
star M150, the agreement between the mass-loss rate given 
by Eq. (fTTj) 7.5 x 10~^° year~^ and the mass-loss rate 
6.4 X lO"^°M0year"i predicted for Z = 0.02 by our hy- 
drodynamic wind models shows a reliability of Eq. 

The winds of ZAMS stars are very weak (with the mass- 
loss rate of the order of 10~^° Mq year~^ for the solar chem- 
ical composition) due to their relatively low luminosity, 
whereas the winds of supergiants may be stronger with the 
mass-loss close to 10~^ Mq year~^. 

For these stars the fitting formula Eq. ([23]) only gives re- 
hable predictions for M > 10^^ M© year , for winds with 
smaller M the error could be of a factor of 3. 

5.3. Stars with 28 000 K < Teff < 40 000 K 

For stars 28 000 K < Teff < 40 000 K, the wind is possible, 
and carbon and nitrogen become more important for the 
driving of the wind, especially for cooler stars. The metal- 
licity needed to drive a wind is about 10"'* — 10"'^ for ZAMS 
stars and may be extremely low (~ 10^^) for most massive 
evolved stars. This is connected with the relatively large 
factor (M/lM0)(i?*/lR0)"^ in Eq. ^ for ZAMS stars 
(these stars are relatively compact), whereas this factor is 
lower for massive evolved stars with large radii. 



5.4. The coolest stars with 10 000 K < Teff < 28 000 K 

The coolest ZAMS stars with T^s < 20 000 K do not have 
any wind because these stars have a too high surface gravity 
to accelerate the wind. On the other hand, evolved massive 
stars may have winds even for metallicities as low as 10~^. 
For Z « 0.01, Pop III supergiant winds may have relatively 
high mass-loss rates of about 10~® M© year"*, comparable 
to present-day hot star winds. The winds of these stars are 
mainly accelerated by singly and doubly ionized carbon and 
nitrogen. Oxygen does not contribute significantly to the 
radiative force in this temperature range as the resonance 
lines of singly and doubly ionized oxygen have higher fre- 
quencies than the hydrogen ionization frequency, and only 
a relatively weak radiative flux is emitted there. 

For the coolest stars, the fitting formula Eq. 
can only be used for higher mass-loss rates M > 
lO-^Mgyear-i. 

6. Consequences of the wind's existence 

6.1. Influence on the stellar evolution 

Our results show that the winds of first stars driven by 
CNO elements are rather weak compared to the winds of 
present-day hot stars. This is caused by two effects, the 
smaller number of CNO lines compared to, e.g., iron ones, 
and by the large surface gravity of first stars. Consequently, 
the estimated total mass lost due to winds during the stel- 
lar evolution is rather small. For example, assuming that 
the final abundance of CNO elements is a solar one, the es- 
timate of the total mass loss of a star with an initial mass 
Mq = 100 Mq is arround IMq, and that of a star with 
Mq = 20 Mq is about 0.1 Mq. Consequently, such stars 
lose only up to a few percent of their mass via line-driven 
winds. 

In reality, the total mass lost might be even lower be- 
cause CNO elements may occur at the stellar surface only 
in the later phases of the stellar evolution and the surface 
metallicity might be lower than the solar one. 

6.2. Metal enrichment of the primordial halos 

The presence of a relatively small amount of metals (with 
Z K, 10^^ — 10^"*) in the primordial haloes m ight inhibit 
subse quent formation of supermassive stars fsee lLoeb et al.l 
|2008| . for a review). Consequently, if the stellar winds of 
massive stars are strong enough, they can change the initial 
mass function in the primordial haloes. 

If in the halo with the baryonic mass of about 10^ Mq a 
first star forms with the mass of 100 Mq, it might lose mass 
arround 1 Mq due to the line-driven winds. Consequently, 
roughly 10~^Mq of freshly synthesized CNO elements oc- 
cur in such a halo, leading to the mass fraction of heavier 
elements of about 10"^, which is too low to change the ini- 
tial mass function. The metallicity of the primordial halo 
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would be enhanced more significantly only after the consid- 
ered star explodes as a supernova. 



7. Discussion 

7.1. The chemical composition 

The mass-loss rate predictions were calculated for the same 
number density of each of the CNO elements. In reality, the 
fraction of these elements is not uniform as assumed here. In 
such a case an approximate estimate of the wind mass-loss 
rate can be obtained as for different temperatures the wind 
is accelerated by different elements. For the hottest stars 
with Toff ^ 40 000 K, the wind is accelerated mainly due 
to oxygen lines (as supported by our models driven purely 
by oxygen), consequently for these stars the mass-loss rate 
depends mainly on the oxygen abundance. For cooler stars 
the wind is accelerated because of both carbon and nitro- 
gen, so it could be roughly assumed that the wind mass-loss 
rate scales either with carbon or nitrogen abundance. 

Stellar winds of present luminuous hot stars are mainly 
accelerated thanks to numerous heavier element lines (es- 
pecially of iron and also o f, e.g . , nickel or coppe r, e.g., 
iPauldrachI [19871: iPuls et al.l[200a IVink et al.ll200l in the 
region where the wind mass-loss rate is determined, i.e. 
below the wind critical point. CNO elements do not have 
many strong lines available to drive a wind as do heav- 
ier elements. In the formalism of line-strength distribution 
functions, CNO elements and iron-peak elements dominate 
in dif ferent parts of t he line-strength distribution func- 
tion (|Puls et all I2000f) . As the result of this, CNO ele- 
ments in present hot stars are important in the case when 
their large abundance (compared to iron-peak elements) en- 
ables their lines to continue beeing optically thick, i.e., in 
the outer wind regions of luminu ous stars or in the weak 
winds of less luminuous st ars fe.g.. lVink et al.|[200ll : lKrtickal 
l2006HKrticka et al.ll2006l ). These trends are reflected in the 
CNO mass-loss rate predictions presented here. Presented 
mass-loss rate predictions of luminuous stars are signifi- 
cantly lower than those derived assuming solar mixture 
of heavier elements (e.g:..IVink et al.ll200H[K]Idritzkil l2002t 



iKrticka fc Kubatll2004lGrafener fc Hamannll2005D . On the 
other hand, the pure CNO mass-loss rate of less luminu- 
ous stars (or stars with very low metallicity) corresponds 
roughly to what is calculated assuming solar mixture of 
elements. 



7.2. Comparison with literature 

Since the wind models driven purely by CNO elements were 
not, to our knowledge, calculated for the studied stars, a 
direct comparison between our results and the results avail- 
able in the literatur e is not po ssible. Pure CNO wind models 
were calculated by lUnglaubl (f200&) . but for different stel- 
lar parameters. However, at a very low wind density, the 
stellar winds are accelerated mostly by CNO lines and the 
contributio n of hea vier elements is ba sically negligible (e.g., 
IVink et aLll2001l : lKrti5ka et al.ll2066l ). Consequently, an in- 
direct comparison between our predictions and the predic- 
tions available in the literature for the lowest metallicities 
is feasible. 

iKudritzkil (|2002D provided mass-loss rate predictions for 
very massive O stars at very low metallicities. The selected 
stellar sample overlaps with ours for the stars with mass 



100 M(7). For the se stars the lowest metallicities for which 
IKudritzkil ()2002f) provides the mass-loss rate predictions 
reasonably agree with our minimum values of the metal- 
licity needed to drive a wind (see Fig. [T]). There is also 
reasonable agreement between our models an d the mod- 
els wi th th e lowest mass-loss rate c alculated by IVink et all 
(I200I and IVink fc de Koted (|2005h . 

7.3. Multicomponent effects and pure metallic wind 

The hydrogen and helium components of the stellar wind 
are accelerated mainly because of the collisions of hydro- 
gen and helium atoms with heavier elements. Consequently, 
stella r winds of hot stars have a m ulticomponent nature 
(e.g., lOwocki fc Pldil2002l : IKrticka 2006). For low-density 
winds the transfer of momentum between radiationally ac- 
celerated metals and hydrogen and helium becomes ineffi- 
cient. In this case the frictional heating or the decoupling 
of wind components may influence the wind structure. 

If the decoupling occurs for velocities that exceed the 
escape speed, then the multico mponent effects do no t in- 
fluence the mass- loss rate (e.g., IVotruba et al.l [20071 ). On 
the other hand, if the decoupling occurs close to the star or 
in the stellar atmosphere, th en the existence of pu lsating 
shells (jPorter fc Skouzalll999f) or pure metallic wind (jBabeli 
|l996 l) is possible. 

Following [Krtickal (|2006D . we calculated the nondimcn- 
sional velocity difference Xhp between a given heavier el- 
ement h and the hydrogen-helium component. For small 
velocity differences Xhp ^ 0.1, the wind is well-coupled 
and can be modelled as a one-component one. For larger 
velocity differences 0.1 < Xhp ^ 1, the frictional heating 
influences the wind temperature, and for large velocity dif- 
ferences Xhp ^ 1, the wind components decouple. The tests 
showed that, for nearly all wind models of ZAMS stars and 
for the low metallicity models of evolved stars, the multi- 
component effects are indeed important. However, in most 
cases the multicomponent effects start to be important in 
the outer wind regions, so these effects do not affect the 
wind mass-loss rate. 

For the cases when homoge neous wind i s not possible, 
a pure metallic wind may exist (|Babe]||l996r ). We postpone 
the study of such stellar winds to a separate paper. 



7.4. Model simplifications 

The radiative force in the models presented here is cal- 
culated neglecting line overlaps and the influence of the 
line transitions on the continuum radiative transfer (see 
IKrticka fc KubafI I2004L for a more detailed discussion of 
the simplification of our code). However, these effects do 
not significantly influence the results here. First, CNO el- 
ements do not have as many lines as iron-peak elements; 
consequently, the line overlaps in CNO winds do not play 
as significant role as in the winds driven also by heavier 
elements. Second, the calculated mass-loss rate is relatively 
low, so the number of strong optically thick lines is rela- 
tively even lower. Finally, even for the winds driven by the 
solar mixture of elements our models are able to predict 
reliable wind parameters. 

Hot star winds display small-scale inhomogeneities 
(clumping ) that may influence t he predicted wind mass-loss 
rate fe.g.. Ide Koter et al.ll2008[ ). As the effect of these in- 
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homogeneities on the wind mass-loss rate is not yet known, 
we neglect this effect. 

7.5. H-He winds and other ways how to lose mass 

A question might appear as to why the stellar wind driven 
purely by hydrogen and helium lines is not possible. This 
question is even more appealing in view of the fact that, 
in some models, the radiative force due to hydrogen may 
contribute a few percent to the total radiative force. To ob- 
tain a wind, in such a case, one would have to decrease the 
wind density. However, as the neutral hydrogen originates 
from the recombination, the decrease in the wind density 
leads to the decrease in the neutral hydrogen fraction and, 
consequently, to the decrease in the radiative force due to 
hydrogen. Consequently, as a result of this coupling the 
winds driven purely by hydrogen are not possible. A simi- 
lar effects also occur for helium. 

We note that there are several other possibilities how 
massive stars can lose their mass, although these possibili- 
ties were not up to now studied in greater detail. Stars with 
the luminosity higher than the correspondi ng Eddington 
limit may lose mass due to porous winds towocki et al 



2004 ) or via 77 Car type of explosions (Smit h fc Owocki 
20061) . Moreover, stars rotating with the critical rota- 



tion rate may lose their mass via equatorial disc (e.g. 



%y lose 
1 120061) . 



iMevnet et al 



8. Conclusions 

We have studied the stellar winds of massive first stars to 
show that pure hydrogen-helium, hot first stars do not have 
any stellar wind. Only stars very close to the Eddingto n 
limit may have very weak pure hydrogen wind (|Paper Ij ). 

As soon as CNO elements are synthesized in the stellar 
core and transported to the surface via the mixing pro- 
cesses, the CNO driven wind may exist for stars that fulfil 
the wind condition Eq. psp . We calculated the models for 
this wind and provided an approximate formula for calcu- 
lating the mass-loss rate as a function of stellar parameters. 

With decreasing Z, both the mass-loss rate and the 
number of optically thick lines decrease. We have shown 
that the condition that at least one line is optically thick 
Eq. is nearly the same as the wind condition, Eq. 
according to which the radiative force should be greater 
than gravity. Consequently, the wind ceases to exist for de- 
creasing Z both because the radiative force is too low and 
because all lines become optically thin. 

We discussed the infiuence of CNO driven stellar winds 
on the stellar evolution and circumstellar environment. As 
the CNO elements are not able to accelerate the wind as 
eflSciently as heavier elements, the wind mass-loss rates 
of CNO winds are relatively low. Consequently, the CNO 
winds do not significantly influence the stellar evolution or 
the circumstellar and interstellar environments. 

We conclude that line-driven winds strong enough to 
influence the evolution of stars far away from the Eddington 
limit did not occur before the first supernova explosions. 
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Table A.l. The wind parameters of studied stars calculated using hydrodynamic models 
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M200-2 


1 


X 


10" 


-4 


7 


.0 


X 


10" 


-UU 


1000 


M200-2 


3 


X 


10" 


-4 


2, 


.0 


X 


10" 


-08 


660 


M200-2 


1 


X 


10" 


-3 


4, 


.2 


X 


10" 


-08 


780 


M200-2 


3 


X 


10" 


-3 


5, 


.5 


X 


10" 


-08 


1490 


M200-2 


1 


X 


10" 


-2 


6 


,4 


X 


10" 


-08 


2650 


M200-2 


2 


X 


10" 


-2 


6, 


.9 


X 


10" 


■08 


3360 


M500-1 


1 


X 


10" 


-4 


8, 


.3 


X 


10" 


-09 


1510 


M500-1 


3 


X 


10" 


-4 


2, 


.7 


X 


10" 


-08 


2440 


M500-1 


1 


X 


10" 


-3 


6, 


.3 


X 


10" 


-08 
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3 


X 


10" 


-3 


9, 


.9 


X 
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-08 
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1 


X 


10" 


-2 


1, 


.3 


X 
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-07 
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10" 


-2 
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.6 


X 
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.6 


X 
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-3 
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.0 


X 
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.9 
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3 
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4, 
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3 
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6, 


.1 
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10" 


-07 
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-3 


1, 


.7 


X 
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-06 
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3 
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10" 


-3 
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X 


10" 


-06 
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1 
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10" 
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3, 
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3, 
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1 
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3, 
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.1 


X 
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-07 
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1 


X 
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.3 
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-07 
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-08 
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10" 


-4 
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.6 


X 
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10" 


-3 


1, 


.3 


X 
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-07 
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-3 


1, 
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-07 
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-2 


1, 


.9 


X 


10" 


-07 
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